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Abstract 

We look at the properties of clusters of order parameter <j) at critical 
points in thermal systems and consider their significance to statistical- 
mechanical ground rules. These properties have been previously ob- 
tained through the saddle-point approximation in a coarse-grained 
partition function. We examine both static and dynamical aspects of 
a single large cluster and indicate that these properties fall outside the 
canonical Boltzmann-Gibbs (BG) scheme. Specifically: 1) The faster 
than exponential growth with cluster size of the space-integrated <f> 
suggests nonextensivity of the BG entropy but extensivity of a q- 
entropy expression. 2) The finding that the time evolution of <p is 
described by the dynamics of an intermittent nonlinear map implies 
an atypical sensitivity to initial conditions compatible with (/-statistics 
and displays an 'aging' scaling property. 3) Both, the approach to crit- 
icality and the infinite-size cluster limit at criticality manifest through 
a crossover from canonical to (/-statistics and we discuss the nonuni- 
form convergence associated to these features. 
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1 Introduction 



An unanticipated - and remarkable - relationship between intermittency 
and critical phenomena has been recently suggested [T] |2j. This develop- 
ment brings together fields of research in nonlinear dynamics and condensed 
matter physics, specifically, the dynamics in the proximity of an incipiently 
chaotic attractor [3] appears associated to the dynamics of fluctuations of an 
equilibrium state with well-known scaling properties 4j. Here we examine 
this connection in some detail with special attention to several unorthodox 
properties, such as, the extensivity of entropy of fractal clusters, the anoma- 
lous - faster than exponential - sensitivity to initial conditions, and the aging 
scaling features in time evolution. 

Our interest is on the local fluctuations of a system undergoing a second 
order phase transition, for example, in the Ising model as the magnetization 
fluctuates and generates magnetic domains on all size scales at its critical 
point. In particular, the object of study is a single cluster of order param- 
eter at criticality. This is described by a coarse-grained free energy or 
effective action, like in the Landau-Ginzburg- Wilson (LGW) continuous spin 
model portrayal of the equilibrium configurations of Ising spins at the critical 
temperature and zero external field. As we shall see, a cluster of radius R 
is an unstable configuration whose amplitude grows in time and eventually 
collapses when an instability is reached. This process has been shown PQ j2] 
to be described by a nonlinear map with tangency and feedback features [3] , 
such that the time evolution of the cluster is given in the nonlinear system 
as a laminar episode of intermittent dynamics. 

The method employed to determine the cluster's order parameter profile 
0(r) assumes the saddle-point approximation of the coarse-grained partition 
function Z, so that 0(r) is its dominant configuration and is determined by 
solving the corresponding Euler-Lagrange equation. The procedure is equiv- 
alent to the density functional approach for stationary states in equilibrium 
nonuniform fluids. The solution found for <p(r) is similar to an instanton in 
field theories , and from its thermal average (evaluated by integrating over 
its amplitude <fto, the remaining degree of freedom after its size R has been 
fixed) interesting properties have been derived. These are the fractal dimen- 
sion of the cluster jE] [7] and the intermittent behavior in its time evolution 
PP [2]. Both types of properties are given in terms of the critical isotherm 
exponent 5. 

As we describe below, the dominance of 0(r) in Z depends on a condition 
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that can be expressed as an inequality between two lengths in space. This 
is ro ^> R, where r is the location of a divergence in the expression for 0(r) 
that decreases as an inverse power of the cluster amplitude 4>o- When r R 
the profile is almost horizontal but for ro ~ R the profile increases from its 
center faster than an exponential. It is this feature that gives the cluster 
some atypical properties that we present and discuss here. These properties 
relate to the dependence of i) the number of cluster configurations on size R, 
and ii) the sensitivity to initial conditions £t of order-parameter evolution on 
time t. 

The above-mentioned properties appear to be at odds with the usual 
Boltzmann-Gibbs (BG) statistics but suggest compatibility [Hj with one of 
its generalizations known as g-statistics jH] [ID]. A condition for these prop- 
erties to arise is criticality but also is the circumstance that phase space has 
only been partially represented by selecting only dominant configurations. 
Hence, the motivation to examine this problem rests on explaining the phys- 
ical and methodological basis under which proposed generalizations of the 
BG statistics may apply. For this reason, we begin in Section 2 by recalling 
the basic framework for the Landau approach, and, in Section 3, we review 
details of the derivation of the dominant configuration. Then, in Section 4 
we consider an estimate for the numbers of configurations that give rise to 
a cluster and check on the extensivity of entropy expressions. In section 5 
we review the connection between the cluster time evolution and the non- 
linear dynamics of intermittency, so that in Section 6 we can discuss again 
the manifestation of g-statistics. In Section 7 we summarize and discuss our 
results. 

2 A Landau approach for single clusters 

We start our account by recalling that the Landau approach is backed by a 
two-stage calculation strategy for the partition function Z of the system. In 
the first, a sum is made over the microscopic configurations that lead to a 
specific form for the order parameter <p(r) (where r is the spatial position) 
to obtain a partial, 0(r)-dependent, result Z^. The second stage consists of 
summing up Z^ over all possible forms 0(r), so that 




(1) 
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where a path integration over the configurations for <p is indicated. The 
partial partition function is then written as 



^ = fi[0]exp(-£[0]/fcT), 



(2) 



where E[cp] is the energy of the system when the order parameter takes the 
fixed form 0(r), and is the number of microscopic configurations for that 
0(r). Taking logarithms to both sides of Eq. (j2J) gives 



where is the so-called LGW free energy and S[<p] = In Q[<p} is an 'entropy' 
term associated to Q[4>] (Boltzmann's constant has been taken as unity). In 
practice the evaluation of 



is considered with ^f[(j)] written as a spatial integral over some function of 
<p(r) that involves a square gradient (perhaps higher-order derivatives) and 
powers of 4>(r). 

At criticality the LGW free energy takes the form 



5 is the critical isotherm exponent and d is the spatial dimension (5 = 5 for 
the d = 3 Ising model with short range interactions). The coefficients a and 
b are dimensionless couplings and <p — A _< V and r = Xr are dimensionless 
quantities, A is the ultraviolet cutoff that fixes the coarse graining scale, <ft is 
the order parameter (e.g. magnetization per unit volume) and r represents a 
spatial position coordinate with dimension of length. Out of criticality 
takes the form \l>[0] = j dr d [l/2 (V0) 2 + 1/2 r o 2 + «o0 4 ] , where r = a t, 
a and u > are constants and t is the reduced temperature distance to the 
critical point. The ordering field has been set to zero [TTj . 

To study the coarse-grained configurations of a single cluster, a subsystem 
of finite size R is considered in the neighborhood of r = and there are no 
restrictions imposed to the value of <j) at its boundaries with the rest of the 
(infinite-sized) physical system. In this case the integration in Eq. is 
performed over the subsystem's volume V. 





(3) 
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3 Cluster properties from dominant configu- 
rations 



We consider a one dimensional system with unspecified range of interactions, 
as in this case the procedure is more transparent and the expressions for 
the relevant quantities more easily derived. The analysis can be carried 
out for higher dimensions with no further significant assumptions and with 
comparable results The LGW free energy reads now 



R 



dx 



dx 



1 5+1 



(5) 



and we adopt the saddle-point approximation - valid for a ^> 1 - to circum- 
vent the nontrivial task of carrying out the path integration in Eq. (|T|). The 
saddle-point configurations are obtained from the Euler-Lagrange equation 



d 2 



dV 



dx 2 



(6) 



where V = — b |0| 5+1 . It is helpful to recall the classical motion analog of this 
problem, a particle in time x under a potential V. Integration of Eq. ([B]) 
yields 



U 



dx 



-b\<f>\ 



s+i 



(7) 



where the constant U is the total particle's energy. Subsequent integration 
of Eq. (JJJ) with U — leads to profiles for critical clusters of the form pQ-[7] 



where 



and 



<fi(x) = A \x — x 



A 



-2/(5-1) 



(8) 



6/2 (S - 1) 



-2/(5-1) 



Xq 



6/2(5-1) O - 



-i 



-(5-l)/2 



(9) 



where x is a system-dependent reference position and cf> = 0(0). The value 
of at the edge of the cluster is 0# = 4>(R) and the cluster free energy is 



R 



2ab / dx 0(x) 



5+1 



(10) 
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This family of solutions give the largest contributions to Z and are the analogs 
of instantons in semiclassical quantum theories, i.e. configurations that min- 
imize the Yang-Mills action 

Similar solutions are obtained for small U ~ pQ [2] where now the 
position of the singularity xq depends also on U. These solutions enter Z 
with a weight exp(— aR \U\) and therefore their relevance diminishes as \U\ 
increases. 



4 Extensivity of critical cluster entropy 

We comment now on the nature of the profile (f>(x) in Eq. (jSJ). This function, 
as well as those solutions for U ~ 0, can be rewritten in the form 

(f)(x) = <f) exp q (kx), (11) 

where exp 9 (x) is the g-exponential function exp (x) = [1 — (q — l)^] -1 /^ 1 ) 

with q = (1 + 5)/2 and k = y/Zb(f)£~ 1)/2 . Because 5 > 1 one has q > 
1 and <p(x) grows faster than an exponential cLS X ^ Xq and diverges at 
xo- It is important to notice that only configurations with R xq have 
a nonvanishing contribution to the path integration in Eq. (^) pQ [2] and 
that these configurations vanish for the infinite cluster size system. There 
are some characteristics of nonuniform convergence in relation to the limits 
R — > oo and xq — > oo, a feature that is significant for our connection with 
(/-statistics. By taking 5 = 1 the system is set out of criticality, then q = 1 
and the profile <p(x) becomes the exponential <fi(x) = <poexp(k x), k Q = ^fa^i. 
(Note that in this case Xq — >• 0). 
The quantity 

$(R) = [ dx(f){x), (12) 
Jo 

or total 'magnetization' of the cluster, is given by 

$(i?) = $ { [ex Pg (kR)] 2 ~ q - l} , R < x , (13) 

where $o =sgn(3 — 5) [20 o / (8 — 3)fc] with q = (l + 8)/2. (We have not shown 
the special case 5 = 3). For 5 = 1, one has $(R) = (pok^ 1 ^[exp (k R) — 1}. The 
rate at which $(-R) grows with R, d$(R) / dR, is necessarily equal to 4>Ri the 
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(14) 



while for 5 = 1 it is d&(R)/dR = <po exp(koR). 

The expressions above may be used to estimate the dependence on clus- 
ter size R of the number of microscopic configurations that make up 
the partial partition function in Eq. (J2J) for the dominant coarse-grained 
4>(x). This dependence may be obtained in a way analogous to that of how is 
obtained the dependence with time of the number of configurations Q for an 
ensemble of trajectories in a one-dimensional dynamical system. Here 'tra- 
jectory positions' are given by the values of in microscopic configurations 
and 'time' is given by the cluster size R. Initially adjacent positions stay ad- 
jacent and Q is almost constant but at later times they spread and Q increases 
rapidly. For chaotic orbits the increment is exponential j3] but for marginally 
chaotic orbits at the tangent bifurcation Q increases as a g-exponential with 
q > 1 [T^j . The ensemble of trajectories is initially contained in the 
interval [0,<f>o] an d at time R they occupy the interval [O,0#], therefore we 
assume 



The results in the following Sections provide a justification for this choice. 
Then, it is significant to note that the Tsallis entropy 

Q l ~ q - 1 

S q = ln q Q = , (16) 

(where lia q y = {y 1 ~ q — — q) is the inverse of exp q (y)), when evaluated 



obtained from S q when q = 1, when evaluated for Q ~ exp(koR) complies 
also with the extensivity property Si ~ R. 

5 Cluster instability and intermittency 

The profile <f)(x) given by Eqs. (jSJ) or (fTTj) describes a fluctuation of the critical 
equilibrium state of the infinite system with average (4>) = 0. In a coarse- 
grained time scale the cluster is expected to evolve by increasing its amplitude 



Q(R) ~ <p Q l d$(R)/dR = O V 



(15) 
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4> and size R because the subsystem studied represents an environment with 
unevenness in the states of the microscopic degrees of freedom (e.g. more 
spins up than down). Increments in 0o for fixed R takes the position xq 
for the singularity closer to R and the almost constant shape 4>(x) ~ <fio for 
xo ^> R is eventually replaced by a faster than exponential shape <fi(x), while, 
as mentioned, the dominance of this configuration in Z decreases accordingly 
and rapidly so. When the divergence is reached at Xq = R the profile (f>(x) 
no longer describes the spatial region where the subsystem is located. But a 
subsequent fluctuation would again be represented by a cluster <fi(x) of the 
same type. From this renewal process we obtain a picture of intermittency. 
A similar situation would occur if R is increased for fixed 0o- 

Indeed, a link was revealed [T] [2| between the fluctuation properties of a 
critical cluster described by Eqs. (JHJ) and (fTTj) and the dynamics of marginally 
chaotic intermittent maps. This connection was demonstrated in different 
ways in Refs. [I] [2] by considering the properties of the thermal average 

($(i?)> = Z- 1 J D\4>] ${R) exp(-* c [0]), (18) 

for fixed R. When xq ^> R the profile is basically flat <f)(x) ~ </>o, the LGW 
free energy is ty c ~ 2abR<j)Q +1 , and the path integral in Eq. ()18|) becomes an 
ordinary integral over < < 0o- One obtains 

($(i?))^^exp(-^ +1 ), (19) 

where u = 2ab(S + l)/{5 + 2)(5 + 3). 

The procedure that resembles the picture given in the beginning of this 
Section is to consider the value of ($) at successive times t = 0,1,..., and 
assume that this quantity changes by a fixed amount fi per unit time, that 
is 

($m) = (^)+/i- (20) 
Making use of Eq. (|19|) one obtains jH] [Zl for small values of <p t the map 

m = e + ^ + z/^ +1 , (21) 

where the shift parameter is e ~ R~ x and the amplitude of the nonlinear 
term is v — ufi. 

Eq. (}2Tj) can be recognized as that describing the intermittency route 
to chaos in the vicinity of a tangent bifurcation The complete form of 
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the map [T] j2] displays a 'superexponentially' decreasing region that takes 
back the iterate close to the origin in approximately one step. Thus the 
parameters of the thermal system determine the dynamics of the map. The 
mean number of iterations in the laminar region was seen to be related to 
the mean magnetization within a critical cluster of radius R. There is a 
corresponding power law dependence of the duration of the laminar region 
on the shift parameter e of the map For e > the (small) Lyapunov 
exponent is simply related to the critical exponent S j2]. 

6 Intermittency and g-statistics 

At the tangent bifurcation, the intermittency route to chaos, the ordinary 
Lyapunov exponent Ai vanishes and the sensitivity to initial conditions & = 
\dxt/dxi n \ (where Xt is the orbit position at time t given the initial position 
x in at time t — 0) is no longer given by the BG law £ t = exp(Ait) but acquires 
either a power or a super-exponential law ^2] |T3j . In this case £ t is given 
by the g-exponential expression, 

& = exp Q (X Q t) = [1 - (Q - l)A Q t]- V(Q ~ 1} , (22) 

containing the entropic index Q and the g-generalized Lyapunov coefficient 
Aq. In the limit Q — > 1 Eq. reduces to the ordinary BG law. See [T2] 
and [T3| for a more rigorous description of the theory and related issues. 

Assisted by the known renormalization group (RG) treatment for the 
tangent bifurcation the formula for £ t has been rigorously derived 
[13] and found to comply with Eq. (|22)1. The tangent bifurcation is usually 
studied by means of the map 

f(x) = e + x + u \x\ z + o{\x\ z ), v > 0, (23) 

with nonlinearity z > 1 in the limit e — > 0. The associated RG fixed-point 
map x' = f*(x) was found to be 

x' = x exp^z/x 2 " 1 ) = x[l -(z- l)vx z ~ l }- l/{z - l \ e = 0, (24) 

as it satisfies f*(f*(x)) = a~ 1 f*(ax) with a = 2 1// ( z_1 ) and has a power-series 
expansion in x that coincides with Eq. ()23|) in the two lowest-order terms. 
(Above x z ~ l = Ixj 2 " 1 sgn(x)). The long time dynamics is readily derived 
from the static solution Eq. (J24|) . one obtains 

Ux m ) = [!-(*- lX 1 ^ 2 '' 2 " 1 ', (25) 
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and so, Q = 2 - z^ 1 and A Q (x in ) = zvx\~ x [12] [HI When Q > 1 the 
left-hand side (x < 0) of the tangent bifurcation map, Eq. (J23J), exhibits a 
weak insensitivity to initial conditions, i.e. power-law convergence of orbits. 
However at the right-hand side (x > 0) of the bifurcation the argument of the 
g-exponential becomes positive and this results in a 'super-strong' sensitivity 
to initial conditions, i.e. a sensitivity that grows faster than exponential [T3] . 
Comparison of Eq. (J23j) with Eq. (|21j) indicates the simple relation z = 5 + 1. 

There is an interesting scaling property displayed by £t in Eq. ()25j) similar 
to the scaling property known as aging in systems close to glass formation. 
This property is observed in two-time functions (e.g. time correlations) for 
which there is no time translation invariance but scaling is observed in terms 
of a time ratio variable t/t w where t w is a 'waiting time' assigned to the 
time interval for preparation or hold of the system before time evolution is 
observed through time t. This property can be seen immediately in £ t if one 
assigns a waiting time t w to the initial position x in as t w = x\~ z . Eq. (|2*5j) 
reads now 

Ut w ) = [l-(z-l)ist/t w ]-^ z - 1 l (26) 

The sensitivity £ t for this critical attractor is dependent on the initial position 
x in or, equivalently, on its waiting time t w , the closer x in is to the point of 
tangency the longer t w but the sensitivity of all trajectories fall on the same 
g-exponential curve when plotted against t/t w . Aging has also been observed 
for the properties of the same map but in a different context |15j . 

7 Summary and discussion 

We have examined the study of clusters at criticality in thermal systems by 
means of the saddle-point approximation in the LGW free energy model PQ j2] 
jH] [I] • The retention of only one coarse-grained configuration leads to cluster 
properties that are physically reasonable but also appear to fall outside the 
limits of validity of the canonical theory. The fractal geometry and the 
intermittent behavior of critical clusters obtained from this method PP [2] (SI 
[7j are both consistent with equivalent properties found for clusters at the 
critical points of the d = 2 Ising and Potts models JH] O n the other 

hand, we found that the entropy expression that provides the property of 
extensivity for our estimate of the number of cluster configurations is not the 
usual BG expression but that of ^-statistics. Likewise, the nonlinear map and 
its corresponding sensitivity to initial conditions linked to the intermittency 
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of clusters do not follow the fully-chaotic trajectories of BG statistics but 
display the features of ^-statistics. 

With regards to the extensivity of entropy, what our assumptions and 
results mean basically is that extensivity of entropy (BG or g-generalized) 
of a cluster and the Kolmogorov- Sinai |3] linear growth with time of entropy 
(BG or g-generalized) of trajectories for the attractor of a nonlinear map are 
equivalent. The crossover from the S q to the Sbg expressions is obtained 
when the system is taken out of criticality, because 5 — > 1 makes q — > 1. 
When at criticality the crossover between g-statistics and BG statistics in 
the dynamical behavior is related to the subsystem's size in the following 
way. We keep in mind that the map shift parameter depends on the domain 
size as e ~ R~ l . So, the time evolution of <fi displays laminar episodes of 
duration < t >~ an d the Lyapunov coefficient in this regime is 

Ai ~ e Within the first laminar episode the dynamical evolution of (f>(x) 
obeys ^-statistics, but for very large times the occurrence of many different 
laminar episodes leads to an increasingly chaotic orbit consistent with the 
small Ai > and BG statistics is recovered. As R increases (R <C Xq always) 
the time duration of the ^-statistical regime increases and in the limit R — > oo 
there is only one infinitely long laminar ^-statistical episode with Ai = and 
with no crossover to BG statistics. On the other hand when R > xq the 
clusters 4>(x) are no longer dominant, for the infinite subsystem R — > oo their 
contribution to Z vanishes and no departure from BG statistics is expected 
to occur [S|. 

In view of the results presented here the departure from BG statistics 
and the applicability of g-statistics is due in part to the presence of the long- 
ranged correlations in space and in time that take place at criticality. These 
correlations give the integrand in the LGW \l/ c a power-law dependence of the 
form l^l^ 1 with 5 > 1 (commonly 5 > 3 as 5 = 3 gives the Gaussian critical 
point) and this in turn determines the g-exponential expression for <fi(x) and 
the properties derived from it. On the other hand, the neglect of all coarse- 
grained configurations other than the most dominant implies that phase space 
has not been properly sampled, and that the ergodic and mixing properties 
characteristic of equilibrium BG statistics are not guaranteed. In this respect 
it is perhaps significant to notice that the most important instance for which 
g-statistics has been to date rigorously proved to hold is the so-called onset of 
chaos in unimodal nonlinear maps JH] This marginally chaotic attractor 
is nonergodic and nonmixing and its dynamical properties exhibit infinitely- 
ranged time correlations. The crossover to BG statistics in this model can be 
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induced via perturbation with noise. Recently j20] , the dynamical properties 
of the noise-perturbed onset of chaos have been proved to be analogous to 
those of supercooled liquids close to vitrification. 
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